Tunable coupling engineering between superconducting resonators: 
from sidebands to effective gauge fields 
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In this work we show that a tunable coupling between microwave resonators can be engineered 
by means of simple Josephson junctions circuits, such as dc- and rf-SQUIDs. We show that by 
controlling the time dependence of the coupling it is possible to switch on and off and modulate 
the cross-talk, boost the interaction towards the ultrastrong regime, as well as to engineer red 
and blue sideband couplings, nonlinear photon hopping and classical gauge fields. We discuss how 
these dynamically tunable superconducting circuits enable key applications in the fields of all optical 
quantum computing, continuous variable quantum information and quantum simulation — all within 
the reach of state of the art in circuit-QED experiments. 

PACS numbers: 



I. INTRODUCTION 

The field of circuit quantum electrodynamics (circuit- 
QED) studies the interaction between artificial atoms 
and artificial photons^'^, both of them implemented with 
the same technology: superconducting circuits cooled to 
millikelvin temperatures. A key feature of these sys- 
tems is that, based on the same microscopic model, both 
the photonic degrees of freedom and the artificial atoms 
have similar energy scales and may interact very strongly. 
Hence, they show effects which are beyond those ex- 
plored in the optical domain. A paradigmatic example is 
the failure of the rotating wave approximation when the 
qubit-photon coupling approaches the qubit and photon 
energies^'^. Aside from the development of qubits and 
the control of their interactions^"^, circuit-QED has re- 
cently started to focus on photons themselves, mostly in 
the context of two different experimental configurations. 
In the first type of setups, cavities are replaced with open 
transmission lines and propagating microwave photons 
that move and interact with localized qubits. This allows 
us to study one-dimensional artificial QED, atom-light 
interaction^, electromagnetically induced transparency^, 
causality quantum metamaterials*"'^"^'^ and to imple- 
ment photodetectors**'^^ and routers*^. The other type 
of setups is based on polariton physics^''' ^^: by coupling 
multiple cavity-qubit systems it is possible to build lat- 
tices on which dressed photons hop and interact, either 
attractively or repulsively, implementing Hubbard type 
models or spin Hamiltonians*^'^°. This gives rise to well 
known models, such as a Tonks-Girardeau gas^*, how- 
ever the architecture based on superconducting cavities 
and Josephson junctions also allows for the exploration of 
new phenomena, such as gauge fields and frustration^^. 

In this work we revisit the architecture of coupled su- 



perconducting cavities, designing a tunable coupling be- 
tween nearest-neighbor resonators. This represents a ma- 
jor breakthrough for this type of systems, because the dy- 
namical tunability of the resonator coupling makes it pos- 
sible to engineer a huge variety of photon-photon interac- 
tions: from red and blue sidebands to gauge fields, pass- 
ing through correlated photon hopping and Kerr non- 
linearities, or simply canceling the usual cross-talk be- 
tween resonators. All these approaches are based on 
simple superconducting circuit elements, such as dc- and 
rf-SQUIDs. 

The setup that we have in mind consists of an array 
of linear resonators, connected through different types of 
Josephson junction (J J) circuits [cf. figure 4]. The fixed 
circuit structure is associated with static, geometry de- 
pendent capacitive and inductive couplings between the 
resonators, while those related to the Josephson junction 
circuits can be tuned by an applied magnetic flux. Most 
notably, the total coupling strength can be reduced and 
even completely suppressed as the coupling due to the 
Josephson junction circuits can have opposite sign. 

Unlike previous proposals that rely on dispersive 
coupling via qubits^'^'^"* or through circulator-like ele- 
ments'^, the coupling scheme proposed here has the po- 
tential to be stronger — the coupling elements are gal- 
vanically coupled, as in Ref. [26] — , it is more robust 
against external perturbations (charge noise) and the ge- 
ometry of the lattice is not tied to the coupling element. 
The coupling elements can be operated in two ways: i) 
with a stationary configuration of magnetic fields that 
determines the associated coupling matrix between os- 
cillator modes or ii) with a periodic multicolor driving 
that allows for engineering arbitrary sideband interac- 
tions, r]ia)h + r]2ab + H.c. between any two resonators, a 
and b, with adjustable coupling strengths 771.2- Moreover, 
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our design achieves tunability in a regime in which the 
couphngs are strong or ultra-strong — larger than the 
corresponding decay rates or comparable to the cavity 
frequencies, respectively — , regardless of other elements 
that may coexist with the cavities, such as qubits or mag- 
netic impurities. 

As potential applications of this work we would like 
to address two fields. The first one has been sketched 
above: by tuning the coupling between different cavi- 
ties it is possible to tune the lattice topology, the cou- 
pling strength and even the phase of the hopping terms in 
polariton arrays. This nicely complements existing pro- 
posals which show how to tune the photon nonlinearity 
by manipulating the qubit inside the cavity^^ and gives 
access to effective gauge fields without relying on frag- 
ile coupling elements^^. The second type of application 
points along the line of quantum information and the ma- 
nipulation of continuous variable states. By means of the 
coupling circuits in this toolbox one may implement any 
nearest neighbor quadratic Hamiltonian with any time 
dependence and geometry, as far as it is embeddable in 
a 2D manyfold. This can be used to implement interest- 
ing states, such as two-dimensional continuous variable 
Gaussian states^'', whose tomography could be supple- 
mented by embedded qubits^^ or moving probes^^. 

The paper is organized as follows. In the first part 
[Sect. II] we will study two superconducting resonators 
that are close together and subject to a mutual inductive 
and capacitive interaction. Using the Lagrangian quan- 
tization, we will show that, both in the weak and strong 
coupling regime, the geometric crosstalk gives rise to a 
constant beam-splitter type interaction. In the second 
part of this work (Sect. Ill) we propose two quantum 
circuits that dynamically tune the inductive coupling be- 
tween the resonators. The first one is a SQUID which 
is galvanically coupled to two resonators. The second 
one uses instead two coupling wires, creating an inter- 
ference device between resonators. We will discuss ana- 
lytically both models, demonstrating that they can tune 
and switch off the overall resonator-resonator coupling. 
In Sec. Ill C we study the validity of our designs under 
realistic experimental conditions, estimating the coupling 
strengths that can be attained in current experiments. In 
Sect. HID we consider the situation of a time-dependent 
resonator coupling. We show that a periodic modulation 
of the coupling makes it possible to engineer sidebands in 
a non-perturbative fashion, controlling the strength and 
phase of both the rotating and counterrotating terms. Fi- 
nally, in Sec. IV we summarize our results and suggest a 
large set of potential applications, ranging from quantum 
information to quantum simulation. 



II. STATIC COUPLING 

In this section we derive the Hamiltonian that rules 
the dynamics of two coupled superconducting strip line 
resonators, and give a general expression for the different 
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FIG. 1: a) Sketch of the geometrical arrangement of the two 
coupled superconducting stripline resonators of length 21. A 
finite interaction is present only in the coupling region of 
length 2lc and is negligible elsewhere, b) Schematics for the 
lumped circuit equivalent. We explicitly draw the mutual ca- 
pacitances (dashed blue lines) and the mutual inductive cou- 
pling (red arrows). The node flux is also indicated. 



coupling constants that arise from the model. Firstly we 
will consider the simplest case of coupling, caused by the 
mutual inductance and mutual capacitance due to the 
spatial proximity of the resonators. Since the coupling is 
time independent and determined by the detailed spatial 
arrangement of the resonators, we refer to it as static 
geometric coupling. The discussion of this interaction is 
done for a particular configuration of parallel resonators, 
but the objective is just to exemplify how this coupling 
manifests itself as a beam-splitter interaction. 

Consider two superconducting stripline resonators of 
length 21, as it is depicted in Fig. la. In this particular 
layout, the coupling occurs mainly within a middle sec- 
tion of length 21 c, where the resonators approach each 
other. Assuming that the crosstalk is given by the mu- 
tual inductance Im and mutual capacitance Cm induced in 
this middle region'^"'^^, we can write down the following 
Lagrangian density: 



(1) 



where both the flux fields, 4>i{x), the capacitance, c, and 
inductance matrices, depend on the position along the 
transmission line 



c = co{x) + c„(.t)(1 - a^), 
I = Iq{x) + lmix)ax , 



(2) 
(3) 



where is the Pauli matrix. A full derivation of (1) 
can be obtained from the lumped circuit equivalent of 
the strip lines [Fig. lb] and it is thoroughly discussed in 
App. A. 
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For the sake of simplicity, we will consider the capaci- 
tance and inductance per unit length of each line, co{x) 
and lo{x), to be constant, and use piecewise constant 
functions for the mutual inductance and capacitance 



a) 







otherwise ' 



\x\ < Ic 

otherwise 



(4) 

We derive a normal mode expansion for the flux 
4>j{x,t) = J2n1j-ni^)''^n{x) in each resonator j — 1,2. 
In what follows we restrict ourselves to the fundamental 
mode of each resonator with frequency ujq and total ca- 
pacitance Cr = j'_iCo{x)dx. Within this subspace and 
mode expansion, the interaction term gives rise to off- 
diagonal terms, as expected from an interaction between 
two cavities, but also diagonal terms that induce a renor- 
malization (dressing) of the oscillator frequencies. This 
dressed resonance frequency is 



ujo\ 



' 1 + C 
1 + 2C 



1 L2 



1 -L2 



(5) 



expressed in terms of two overlap integrals 



uo{x)'^dx, A2 



[da;UQ{x)f dx, (6) 



where C = Cm^i/Cr, L = Im/lo, and v — uj'^Crlo/ ^2 is 
a geometric factor. 

We finally proceed with the quantization of this 
model, introducing the oscillator length oq ~ 
^Jh{\ + C)/CrU!{l + 2C). We express the phase space 
operators in terms of the Fock operators, qj — aQ^aj + 
aj)/\/2 and pj = ih{aj — aj)/\/2aQ. This leads to, 

H = huj a^-a^ — hgc{a\ — a^){a\ — (7) 

- hgi{a\ +ai)(4 +02)- 

The coupling constants gc and gi account for the static 
capacitive and inductive contributions to the coupling, 
respectively, 
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The usual limits in quantum optics correspond to the 
weak coupling and strong coupling regimes. In both of 
them gc^i/bjQ ^ 1, so that the frequency renormalization 
becomes negligible (provided that C,L <C 1). We can 
then invoke the Rotating Wave Approximation (RWA) 
and transform (7) to the beam splitter model. 



K9t+ 9c){a\a2+a\a-^). (10) 
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FIG. 2: (color online) Architectures leading to a tunable mi- 
crowave beam splitter, a) A dc SQUID (superconducting loop 
interescted by two Josephson junctions) mediates the coupling 
between two stripline resonators. The established pointlike 
contact between the resonators takes the coupling to the ul- 
trastrong domain, b) A superconducting ring intersected by 
a Josephson junction now partially shares its branches with 
the cavities, improving the switching capability. 



Note how the resulting Hamiltonian can be interpreted 
as an exchange or hopping of excitations between modes, 
similar to optical lattice and tight-binding models. 

This type of static geometric coupling is implicit in 
the experimental configurations of coupled cavity mod- 
els^'^'i^, though previous designs have inclined to con- 
sider a capacitive coupling taking place at the electric 
field nodes (current anti-nodes) of the resonator-'^^'^^-^^^ 
sometimes enhanced by an additional J J circuit^^'-^s 



i=i,2 



III. TUNABLE COUPLING 



In this section we study alternative mechanisms for 
coupling two or more linear resonators. On the one hand 
we aim at a larger coupling strength, and on the other 
hand we wish to achieve real-time tunability of the cou- 
plings. For both goals it will be advantageous to rely on 
inductive rather than capacitive coupling. First of all, the 
inductive coupling realized by JJs and loops intersected 
by JJs (e.g. SQUIDs) can be tuned by an applied mag- 
netic field varying the magnetic flux threading the JJs or 
loops. Second, and equally important, inductive inter- 
actions can be enhanced, proflting both from the kinetic 
inductance of thin superconducting films and from em- 
bedded junctions working in the linear regime'^''. Based 
on the two previous ideas, we envision the two coupling 
elements sketched in Fig. 2. We will study analytically 
these designs, deriving expressions for the effective inter- 
actions and coupling strengths. 
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A. SQUID as a coupler 

Given the large inductance provided by the Josephson 
junctions, one could naively think of connecting both 
cavities with a superconducting wire interrupted by a 
Josephson junction. In doing so one would achieve a 
static ultrastrong coupling. However, for tuning the 
Josephson inductance of a single junction we have to 
generate a magnetic flux of the order of a flux quan- 
tum threading the junction area. Due to the small junc- 
tion area unpractically large magnetic fields would be 
required. Fortunately, we can design a much better tun- 
ability by using a superconducting quantum interference 
device (SQUID) configuration, as depicted in Figure 2a. 
Since the SQUID loop area is much larger than the junc- 
tion area, much smaller control fields are required. Note 
that even though the use of dc-SQUIDs as tunable cou- 
plers was already suggested for flux qubits'^"', our setup 
works with continuous variables and has subtle differ- 
ences that need to be discussed below. 

A short line with the SQUID represents a small con- 
tribution to the original Lagrangian density (1), that is 
£t = £ + £sQUiD with 

/:SQUID = E%^<^'.,'c + i^.,.COs(^), (11) 

where, 4>j^k represent the flux differences along the junc- 
tions k — 1,2. We use fluxoid quantization along the 
SQUID loop, -I- <j)j,2 + — "'J'o, to express the La- 
grangian in terms of the variables (j)± — 1(07,1 ± ^7,2) 
and the total flux enclosed by the loop, $0. For simplic- 
ity, we assume $0 ~ $ext, that is, we are neglecting the 
additional flux generated by the circulating loop current. 
This is equivalent to restricting our discussion to screen- 
ing parameters (3^ = 2nLIc/^o <C 1 as discussed in more 
detail below. Here, L is the loop inductance and Ic the 
critical current of the Josephson junctions. If the SQUID 
is symmetric, Cj.i = Cj,2 and Ej^i — Ej^2, the coupling 
becomes i?effC0s (7r0_/<I>o), with an effective Josephson 
coupling energy that depends on the flux threading the 
SQUID loop, 

£;eff = 2i;jcos(7r$o/*o). (12) 

The voltage-phase relation ^_ = <j)i{x) — <j)2{x) allows us 
to express </>_ in terms of the voltages at the edges of 
the connecting wire. In the linear limit of small fluxes, 
i.e. small photon number [cf. App. B], we can write a 
quadratic coupling between fields 

-CsQUID ^ ^ (^1 - ^2) - ^^^f°^ (<^i - '^'2)^, (13) 

which by means of the normal mode decomposition 
adopts the form of Eq. (7). The tunability of this model 
relies on the fact that the Josephson energy Egg is flux- 
dependent. This implies that 

47r2 

9^ = 5f + -ttEj cos (7r$o/<i>o) (M) 



in Eq. (10) can be changed in magnitude and sign. For an 
appropiate value of the external flux (close to $0 = $0 /2 
if \gc\ <C \gi\) we can fully deactivate the coupling. 

We finally notice that our setup is robust against small 
differences in the two junction energies, Eji 2 — Ej{l ± 
e). In this case one can still expand (f)i^2 — ^ 4'-i 
linearizing around 0_ ~ to obtain 

£sQUiD±2££;jsin(7r$o/^o)^^- (15) 

<l>o 

The linear term in this equation amounts just to a dis- 
placement of the oscillators and does not add up to the 
total coupling, preserving the tunability of the setup. We 
will use this idea in the following setup. 



B. Superconducting ring coupler 

The second design is shown in Fig. 2b. It consists of 
a superconducting ring interrupted by a single Joseph- 
son junction. Since both resonators share a branch of 
the loop that couples them, the Lagrangian acquires new 
contributions with the kinetic inductance of the super- 
conductor'^'^. This kinetic coupling can be very strong, 
while still retaining the switching capability due to the 
fluxoid quantization inside the loop, similar to previous 
designs for superconducting qubits^^. 

Our derivation is based on two non-essential con- 
straints. The first one is that the loop is small enough 
to neglect its self inductance (/3l ^ 1)- The second one 
is that the wire without junction touches the resonators 
at points where the flux of the coupled modes is zero''^ 
'"o(2;i) — 0. Under these circumstances the coupling term 
reads 

Cjj^^^^ + Ejcos{^-^y (16) 

The fluxoid quantization inside the loop, A(/)2 — A</)i -|- 
(f),j — —$01 allows us to get rid of the flux variable 0j 
and rewrite the coupling in terms of the branch fiuxes 

A0i,2, 

Cjj = ^(a0i-A<^2)' (17) 

At this point we will repeat the linearization of the 
cosine, much like in (13). However, now the Taylor ex- 
pansion will depend on the external flux, 4>f), produc- 
ing linear and quadratic contributions of different mag- 
nitude. We start with the normal mode decomposition 
of the branch fluxes and restrict ourselves to the lowest 
energy modes 

A0J = (j)j{x2) - (t)j{xi) ~ q[^^dxUn{x)\a:^/^X. (18) 



5 



Introducing these terms in the interaction produces the 
quadratic Lagrangian for the fundamental modes, 



2 ^ 

7j(9i - 92) 



ajq,-Pjqj, 



(19) 



ajqiq2 + (3jqiq2- 



The expressions for all coefficients can be computed from 
first principles 



Ej^d^uoix = Ol^Ax^ cos ( 



27r^ , ,^ /27r$o 
7j = -Bj — dj^uoia; = OjAxsm — — 



(20) 
(21) 

(22) 



Out of these terms, 77 is a linear displacement of the 
cavity eigenmodes and does not transfer energy. The 
capacitive and inductive terms, aj and are the only 
ones that contribute to the inter-cavity coupling, gc and 
f/i, and to the frequency renormalization. More precisely, 
we obtain the model (7) with mode frequency 



UJ 




and coupling strengths 



9r 



9c = 9c 



^static 



(23) 



(24) 



In general we will find that for a junction that works in 
the flux regime the term /3j dominates all other contribu- 
tions. But even without this assumption, it is true that 
while gc is fixed, the value of gi depends entirely on /3j 
and can be changed in magnitude and sign, either en- 
hancing the strength of the beam-splitter coupling (10), 
or switching it off entirely for a value of $f) ~ $0/4. 

While the coupling strength grows with the loop size. 
Ax, we can not make it arbitrarily large because then 
we are no longer allowed to neglect the additional fiux 
(pL caused by the circulating loop current due to the in- 
creasing value L of the self- inductance of the loop'^^. In 
this case the total flux threading the loop is given by the 
sum of the external flux "I>ext and the flux However, 
as explained in Ref. [35] chapter 8.4, provided that 



(25) 



we can ensure that the $f) versus $oxt dependence is 
single-valued allowing us full tunability of the coupling. 
This condition means that the maximum loop current 
Ic cannot generate more than a single flux quantum. It 
restricts ourselves to loop sizes of around the 5% of the 
resonator length. We will now study various methods 
to increase the coupling strength while preserving the 
condition above. 
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FIG. 3: Coupling strength for differerent Niobium transmis- 
sion lines, as a function of the loop size. While homogeneous 
resonators (blue) hardly reach the strong coupling regime, in- 
homogeneous ones (red-dashed) do. The coupling strength 
can be further increased with a Josephson junction interrupt- 
ing the center conductor, as shown for the homogeneous case 
(dotted) and the inhomogeneous one (dot-dashed). We have 
considered for each resonator a frequency cjo/Stt — 6.65 GHz. 



C. Estimation of the coupling strength 

We are interested in an upper bound for the coupling 
strength "g" . More precisely we would like to access both 
the strong and ultrastrong coupling regimes. Strong cou- 
pling means that it is possible to observe Rabi oscillations 
between both cavities because the coupling g is larger 
than the resonator decay rate, k. On the other hand, 
ultrastrong coupling occurs when the RWA fails, which 
in this case implies that the number of photons in the 
ground state, that is proportional to g/uj, approaches 
one. 

Looking at the flrst proposal [cf. Fig. 2a and Sect. 
Ill A] we note that the maximum coupling is reached for 
an external flux = n<i>o threading the SQUID loop, 
and thus yielding 



4e^ Ej 
2C' h?uj 



\u,ix)\' 



-Z\ui{x)\' 



(26) 



where Ic is the critical current of the junction, Z is the 
resonator impedance, and the eigenmode uo{x) satisfles 
< |ui(a:^)| < V^- To preserve the power fleld expansion, 
we suppose the SQUID to be built at a position such 
that |ui(a:)| < 0.1 [see Sect. HIE]. Under this condition, 
and using a critical current /j, ~ 5 x lO^^A, together 
with Z = 50 f2, it would be possible to reach a coupling 
strength up to ~ 1.2 GHz. 

On the other hand, for the second proposal [cf. Fig. 2b 
and Sect. IIIB ] the coupling (24) in the loop becomes 



(27) 



with 0x4' = ^^^^/2cf~i^' ^'^^ ^ homogeneous resonator 
(see (Al)), we can straightforwardly assess the slope of 
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ui{x), finding an exact expression for g: 
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2 ^J^c 
Wo 



Ax 

~2r 



(28) 



where ujj — Ej/h, ujc — Ec/h — {2e)'^/2hCr is the 
characteristic charging frequency of the resonator, and 
loq the first-mode frequency. Using available values for 
Nb striplines and junction parameters, we find that the 
homogeneous resonator remains in the weak coupling 
regime, as we envision before. For a loop size Ax of 
1% of the resonator length, we obtain .g ~ 2 MHz which 
represents 0.03% of the resonator frequency lj [see Fig. 3]. 

Adding a constriction to the central part of the res- 
onator increases the field slope and thus the coupling. 
To this end, for a suitable Nb inhomogeneous transmis- 
sion line resonator'^^, this enhances the coupling up to 
g ~ 100 MHz, or equivalently to 1.8% of the resonator 
frequency [Fig. 3 red-dashed]. 

Finally, the coupling can be further enhanced by in- 
terrupting the transmission line with a Josephson junc- 
tion. Due to the presence of the junction, the flux eigen- 
mode presents a constant phase slip A(/)o at a; = 0, 
which depends on the Josepshon coupling energy of the 
junction'^^''^^. This additional phase jump enhances the 
coupling as follows 



^ , ' 27r dtbix) ^ 



A0O 



(29) 



Optimal parameters for the junction attached to the res- 
onator (i^jres — 7i?jioop) could lead to extremely large 
couplings of around g ^ 600 MHz (9% of uj). 



D. Sidebands 

So far we have discussed the possibility of tuning the 
coupling strength between two resonators, constructing a 
classical switch that allows us to control the exchange of 
photons. In this section we discuss a second type of tun- 
ability, which consists of engineering an arbitrary linear 
coupling type between two resonators: 



gie"t'^a^b 



92e 



"f'^ab + R.c. 



(30) 



represented by the Fock operators a and b. This would 
enlarge the applicability of our setup, extending it to the 
realization of almost any quadratic model with nearest 
neighbor interactions. 

In order to demonstrate that this is possible we start 
our discussion by noting that both the dc-SQUID and the 
ring coupler provide us with a flux-dependent coupling 

H = hujaa^a + hujbbh + g{<S>o){a^ + a){b^ +b). (31) 

If we now engineer the two resonators to have very differ- 
ent frequencies, Wq and ujb, a static coupling \g\ ^ uJa.b 



will be effectively suppressed, giving rise to a small dis- 
persive term 

H ^hujaa^a + nwi,b'^b+ -a^ab'^b. (32) 

However, if we allow for a two-tone driving of the cou- 
pling 

g{t) ^ g[^fj{t)] = 51 cos[(tJb - Wa)i + 0i] + (33) 

+ 52COs[(i:Ja +Wf,)i + 02], 

then this driving effectively activates the rotating and 
counter-rotating terms, with the phases given above. To 
show this we switch to an interaction picture with respect 
to the two harmonic oscillators 

Hi = g{t) (a^&e*('^«-""')* + abe-'^"-+'^'''>' + H.c.) . (34) 

The oscillating terms in Eq. (33) will precisely cancel 
the ones in the previous time-dependent Hamiltonian, 
leaving back some other non-resonant terms which only 
act in higher-order perturbation theory, 0{g'^/uj). The 
result should be the desired combination of sidebands 



Heff = gi{a)be"t'^ + abe"''^ + H.c.) + 



(35) 



It is worth mentioning that the previous sideband engi- 
neering is not perturbative: while we still need to impose 
that |gi,2| ^ jwf, — Waj, the resulting coupling is larger 
than the dispersive term. Another very important fea- 
ture is that this method allows us to control the phase of 
the rotating and counterrotating terms, for this is related 
to the phase of the two-tone driving. As we discuss be- 
low, this is a very important property, as it allows us to 
implement effective gauge flelds that control the hopping 
of photons between resonators. Moreover, we achieve this 
effect by a simple driving of a standard SQUID, without 
the need of time-reversal breaking circuits which might 
be very sensitive to other noise sources^'^. 

Finally, even though the realization of the time de- 
pendence (33) might seem very complicated, in practice 
we do not need to tune the flux in a very complicated 
manner. A simple driving of ^oi^) — $ + cos{ujt), 
when introduced in the sinusoidal coupling (14), g ~ 
cos(27r$5/$o) produces, via the Jacobi- Anger expansion 

g{t) ~ cos(27r-I>/$o)^o(<5$)+sin(27r$/$o)^i(<5"I')cos(a;i)-f 

(36) 

in terms of the Bessel functions Jq and Ji. This series 
contains the basic driving plus higher harmonics which 
will be spectrally suppressed in the coupling term. Alter- 
natively, a suitable dependence for can be engineered 
with around 0.1 ns resolution using appropriate signal 
generators. Again, out of this signal only the resonant 
terms, with frequencies around uJa i Wfc will contribute 
to the coupling. Discretization errors in the signal, and 
higher harmonics, will be averaged out. 



E. Nonlinear photon hopping 

So far we have worked with the Josephson junctions in 
the hnear regime, neglecting higher order terms, which 
are of the order ^Ej{2Tr(l>/^o)'^. This approximation is 
vahd only when the argument of the trigonometric func- 
tions, 27r/(/)/$Oj is small, a condition which can be re- 
casted as a restriction on the number of photons that 
can populate the resonator. Roughly, for the SQUID we 
have the condition 

where Z is the impedance of the resonators, n is the 
average number of photons and u{x) is the mode wave- 
function at the coupling points. Using, in the same way 
as above, the value u{x) — 0.1 restricts the number of 
photons to be n < 1000 [see App. Bl], which docs not 
represent a restriction for the few photon applications 
that we envision. 

The question now is what happens when we do not 
neglect the nonlinear terms. In this case we have the po- 
tential to introduce new interactions between resonators, 
which are now of higher order and include on-site non- 
linearities, nf, nearest-neighbor attractive or repulsive 
interactions, niUj, photon-pair hopping, al^a^, etc [cf. 
App B 2] . Out of these terms some are already strongly 
suppressed because of being off- resonant; this is the case 
for interactions with odd powers, such as a] aj. The 
Kerr nonlinearities will always be present and give rise to 
extended Bose-Hubbard physics. Finally, the correlated 
hopping terms, aJ^Oj, can be resonantly enhanced using 
the same technique that we employ for the sidebands: 
introducing a frequency mismatch between neighboring 
cavities and driving with exactly the frequency which is 
needed to select this process, 2{uji — ujj). With all these 
tools we envision the possibility of engineering very inter- 
esting models, such as a condensate of pairs of photons"^^, 
which are very hard to engineer in other systems. 

IV. APPLICATIONS AND DISCUSSION 

Summing up, in this work we have studied two different 
ways to engineer the coupling between superconducting 
resonators: one is geometric and static in nature, while 
the other relies on nonlinear coupling circuits and can be 
easily tuned in and out of the strong coupling regime. 
Both elements together form a powerful toolbox for im- 
plementing almost arbitrary models consisting of a low- 
dimensional (from ID to 2D) array of resonators with 
tunable nearest neighbor interactions, as in the model 
sketched in Fig. 4. Let us now discuss some of the po- 
tential applications of such circuits. 

a. Traditional Quantum Optics. The implementa- 
tion of tunable sidebands in coupled resonators opens the 
door to many well-known processes from quantum optics. 
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FIG. 4: a) One dimensional array of cavities coupled by means 
of a superconducting ring coupler, b) Two-dimensional lat- 
tice of circular resonators, coupled by dc-SQUIDs. Both lat- 
tices are bipartite. Using different resonator frequencies for 
each sublattice, uua ij^b, we can use the techniques from 
Sect. HID to engineer any sideband interaction between the 
arrays. 



Some of them are the squeezing of different modes via 
those sidebands, frequency conversion of photons as they 
are transferred between cavities, parametric generation 
of photons via a^b^ + ab terms, entanglement produc- 
tion at high temperatures^®, etc ... The beam-splitter 
Hamiltonian is also the cornerstone of all-optical quan- 
tum information processing, as suggested in Ref. 25 for 
a different circuit-QED architecture. 

b. Harmonic models. The most immediate applica- 
tion of our design would be to implement arbitrary quasi- 
local and quadratic Hamiltonians, with the aims of study- 
ing the dynamical or statical properties of many-body 
Gaussian states. This includes a variety of studies, such 
as the static correlations in the model'^^ and their rela- 
tion to the underlying entanglement, dynamical quantum 
phase transitions from trivial to critical phases, the study 
of propagation of correlations in non-equilibrium models 
and their relation to Lieb-Robinson bounds'"'"^^ , etc. In 
this context, the tunability of the coupling plays two dif- 
ferent roles. On the one hand it allows us to change the 
parameters of the Hamiltonian in an abrupt or smooth 
way, for instance to study a dynamical quantum phase 
transition or a quench. On the other hand and equally 
important, by switching off all couplings we can freeze 
the quantum state of the oscillators, giving us time to 
measure the properties of the system, either with differ- 
ent measurement qubits or through a movable probe^^. 

c. Anharmonicity & thermalization. The interest of 
the harmonic problems lays in their simplicity and the 
possibility of obtaining analytical and numerical results 
for different geometries and sizes. However, as soon as we 
introduce a small nonlinearity in our system, we can say 



8 



very little about their dynamical and static properties 
and many of the simulations which we mentioned in the 
previous paragraph become open problems. In particu- 
lar, one very simple problem which deserves being stud- 
ied is that of thermalization. The basic idea is to replace 
the linear resonators in Fig. 4 with resonators that host 
a tunable and weak nonlinearity in the form of a SQUID 
(similar to Ref. 43 but outside the linear regime). One 
would then prepare the ground state of the cavities with 
a value of the coupling, and then abruptly quench this 
coupling to a different (larger or smaller) value in which 
the prepared state is not a ground state. Throughout 
this process it will be possible to track the relaxation of 
the oscillator chain or lattice, studying how its behavior 
is modified by the presence of the nonlinearity. 

d. Coupled cavities & gauge fields Along the lines 
of anharmonic systems, another interesting problem is 
the study of coupled cavities or Jaynes-Cummings lat- 
tices^^'^^. The setup would be that of Fig. 4, but with 
one qubit attached to each resonator. The coupled 
qubit-resonator system behaves as a highly nonlinear el- 
ement, implementing a quasiparticle known as "polari- 
ton" , which may hop from resonator to resonator through 
our tunable coupling elements. This can be roughly for- 
mulated as a Bose-Hubbard Hamiltonian 



J7(a|ai), 



(38) 
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Appendix A: LUMPED ELEMENT MODEL OF 
TWO COUPLED RESONATORS 

Here we derive the density Lagrangian (1) of Sect. II 
from the quantum network theory perspective. 

The appendix is divided into three parts: in the first 
one we review the quantization of a single microstrip res- 
onator. In the second one, we consider the equivalent 
circuit of the coupled strip lines shown in Fig. la in its 
lumped element model [see Fig. lb]. The Kirchhoff equa- 
tions derived here will give rise to the Lagrangian (1) in 
the continuum limit 



with a very nonlinear on-site interaction U and a hop- 
ping tij which, unlike previous proposals^'', is now dy- 
namically tunable. This allows us to explore the quan- 
tum phase transitions from weak interactions U <Si \t\ to 
hard-core particles ?7 ^ |t| simply by reducing the hop- 
ping instead of arbitrarily boosting the qubit-resonator 
interaction — something which might be more challeng- 
ing from the theoretical and experimental point of view. 

In addition to the usual Mott-superfluid phase transi- 
tion, we now have control over the phase of the hopping, 
tij — \t\ expiOij. The procedure, as described in the pre- 
vious section, consists of engineering two coupled sub- 
lattices [A and B in Fig. 4b] of resonators with different 
frequencies, loa wb. Applying a multitone driving on 
the bonds that connect both sublattices, we can create 
an array of phases 9ij which have a nontrivial flux around 
each plaquette. This will allow us to probe the integer 
quantum Hall physics with polaritons, without the use of 
circulators^^. 

In summary, we have shown that in circuit-QED, tun- 
able coupling between resonators can be implemented via 
simple Josephson circuits. We have developed this initial 
idea into a profound theoretical basis for exciting multi- 
resonator experiments ranging from arbitrary sideband 
interactions to setups scalable towards the many-body 
regime. On the theoretical side, our results lend them- 
selves to be expanded to advanced scenarios, such as the 
relation between our circuit models and Josephson junc- 
tion arrays, the influence of decoherence, or even for the 
design of models with tunably dissipative elements. 



1. Single oscillator description 

Here, we detail the description for the single resonator 
case. The transmission line held equations are obtained 
from their lumped circuit equivalent. Neglecting losses 
it can be described as a series of LC circuits''^. In the 
continuum limit, the resulting held equations can be ob- 
tained from the Lagrangian: 



Co 



dx 



-I 



co{x)^{x,tf -h{x)-^d^(l){x,tf , (Al) 



where Cq{x) and ^0(2^) ^-re the capacitance and induc- 
tance per unit of length, respectively; otherwise, ipix) — 
(<i>o/27r)<p(a;) is the magnetic flux variable with $0 = 
h/2e the magnetic flux quantum and ^pix) the phase of 
the macroscopic wavefunction describing the supercon- 
ductor. The stationary modes are found by solving the 
eigenvalues and eigenvectors for the equation of motion 
(the Euler-Lagrange equations) 

d,[h{x)-^d^c^{x,t)\ =co{x)d^<f>ix,t), (A2) 

which is nothing but the wave-equation in one dimen- 
sion. The solution to this equation is expanded in normal 
modes and time dependent amplitudes. 



(2;,*) ^^qn{t)Un{x), 



(A3) 



such that g„ = —LUnQm with a;„ the resonator frequen- 
cies. Therefore the eigenstates m„ satisfy the differential 
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equation i9a; [^0 (a;) ^dxUn{x)] = —uJnCo{x)un{x). The w,, 
satisfy the orthogonal relation: 



dx. 



Co,jix)Um,]{x)Unj{x)dx = CrSr, 



(A4) 



Hence, the Lagrangian (A6) ends up in Eq. (1) in the 
main text that we rewrite here for completeness, 



with Cr — J_iCo{x)dx the total capacitance of the res- 
onator. 

For homogeneous resonators, Iq and Cq are constant, 
and we obtain the well known case of equispaced eigen- 
frequencies = (2n — l)/2/-\/ZoCo with 21 the length of 
the superconducting resonator and m„ = \/2sin((27T, — 
l)7ra;/2/). 



C 



1 



E 

i=i,2 



dx 



Cij(j)i{x)4>j{x) - l^ ^dx(j)i{x)dx4>j{x) 



1. Two coupled oscillators 

Consider the lumped element model depicted in Fig- 
ure lb. The discrete modes of the electromagnetic field 
inside the strips are described as arrays of LC oscillators, 
together with the mutual inductances and capacitances, 
representing the crosstalk. By applying the current con- 
servation law at each node of the circuit, we obtain a set 
of dynamical equations for the node fluxes (^j,n'- 

Axco{n)(f)i^n + Axc,„(n)(^i,„ - 02,™) = (A5) 



(A9) 

The fluxes, (t>i{x), are thus coupled by the capacitance 
c{x) and inductance l{x) matrices given in the main 
text (2). The diagonal terms of these matrices represent 
the single resonator Lagrangian Cq derived in the pre- 
vious section, that depends on lo{x) and co{x). On the 
other hand, the off-diagonal contributions represent the 
coupling Lagrangian Ci, described by the parameters 



and C'ij — — 



3. Generalization to more oscillators 

We now show that the quantum description of two cou- 
pled resonators presented above can be generalized to the 
case of TV coupled resonators. We therefore extend the 
sum in (A9) to N 



Ax lo{ny - lm{n)^ Ax Zo(n)2 - -f 1)^ 

and the equivalent equation for the second resonator. 
Above Co (n) , Iq (n) are the capacitance and inductance 
per unit length respectively, while c,n (n) , Im (n) repre- 
sent the mutual capacitance and mutual inductance co- 
efficients. Notice that, in general, all these parameters 
are position dependent. 

The former equations of motion are nothing but the 
Euler-Lagrange equations associated to the following La- 
grangian 



1 



N 

E 



dx 



Cij4)i{x)4ij{x) - lj^j^dx(j)i{x)dxcj)j{x) 



(AlO) 

where the c{x) and l{x) are now given by x iV matri- 
ces with self capacities (self inductances) on the diago- 
nal and the mutual capacities (self inductances) on the 
off-diagonal. Following the same procedure as above we 
restrict ourself to the fundamental modes, split off the 
single resonator Lagrangians and write the interaction 
part as 



L = T-V, 



(A6) 



-^1 = ^E ((^mAiq. 



with 



N 
N-1 



A „2 



(All) 



Ax ■ 



= ^ E + Cm(?l)(0i,n - ^3 + l,n f, (A7) 



( -c,nAiq.iqi+i 



,2 _ ,2 ^2qiqi+l 



-E 

Ax ^ 

^0 (^)'^ 



^=2A;,l.,g37:t(^(^."-^.«-)^ (A8) 



only taking into account nearest neighbor interaction. 
The Hamiltonian can finally be written as 



N 



JV-1 



We can now take the continuum limit Ax — > 0, which 
implies: 



-l) 



Ax 



H/h - ^c^a]a, - - «,)(«]+! - 

AT-l 

with w, (7c a-nd gi identical to the two resonator case. 
In particular from the resulting total Hamiltonian H the 
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normal frequencies can be found, giving: 



1 



W_ = UJQ 1 



1 



L 



UJ+ = Wo W 1 - 



1 + 2C V ^^(1--^) 
' L 



(A13) 



1/(1 + 



We finally notice that by making Ic ^ I { i' I), i.e. two 
straight parallel resonators, the formulas for the normal 
modes match the case of two coupled LC circuits. 



Appendix B: Linear and non-linear couplings 

The motivation in this appendix is twofold. On the 
one hand we estimate the validity of the linear approx- 
imation; on the other hand, we explicitly compute the 
first non-linear corrections to the coupling. 

We first expand the cosine in £squid (H), 



-cos 



2tt 



= -1- 



1 /27r 



We recall that: 

= uq\/ hZ {a[_ + a_) , and a_ 

Assuming Z — bOfl and defining 



1 /27r 
24 



(Bl) 



V2 



(ai-as) (B2) 



f EE — 2hZ = u^ — 
we can write for the expansion: 



(B3) 



-cos 



^0 U-l+f(al+a-r- 



96 



(al+a^)^ (B4) 



1. Linear regime 

The linear approximation is justified when the second 
order terms in Eq. (Bl), or equivalently the average value 



and fluctuations of the flux in Eq. (B4), are small. How 
does this relate to actual experiments? In order to de- 
termine a condition based on the number of photons we 
study the fluctuations (</>?_), which are related to the ex- 
pectation value 



((al + a_)2) 2{{a\ai) + (4a2>) = 4n 



(B5) 



with n the number of photons. Using the previous series 
we conclude that linearization is strictly justified when- 
ever n <^ 10/mq, where uq is the value of the mode wave- 
function. For the ring coupler layout, the same reasoning 
follows by replacing ua — ?> OxUqAx in (B2) and (B3). 
2. Non-linear hopping terms 



With the help of Pathak's results^**, we compute, 



{a^+alf = at +A{alfa_ +6{al)^a'i +6a'' a + S + R.c. 

(B6) 

Notice that a_ = TJ^"^ ~ '^^^ (^^)]' above 

ends up in a long expression hard to deal with. To make 
it simpler, we assume that sidebands will select or a^a 
and we resort to a RWA argument to write: 



(a_+al)* Gialfal + Ua"^ a + 6{alf + 6al+3. (B7) 



Looking at each term we rewrite in terms of the local 
bosonic operators, ai and 02: 



4{alfa^_ = {a\)^al + {alfaj + 2a|4aia2 (B8) 
-2a|4 (a? + 4) + H.c. 
2{al)^ = (4)2 + (4)2 -244 (B9) 
2a^a_ = a|ai + 0302 — a|a2 — a2ai. (BIO) 
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